Magnetic field data, including square maps used for figures and transects used for peak field analysis.
Outputs are organized into three groups:
1. BMaps: Final magnetic field maps
2. BMaps_nolr: Magnetic field maps neglecting late remagnetization
3. BMaps_no2k: Magnetic field maps without the top layer of magnetization (set by the resolution of the magnetization models; corresponds to top 1 km for basins <1000 km in diameter, top 2 km for 1000 km diameter basins, top 4 km for 2200 km diameter basins)
4. RevRates: Mean reversal rates

Terminal folder names are typically of the format [polarity vector]_[imaging grid used]; for example, the folder name 010_ig3 assumes a polarity vector oriented in the (X,Y,Z) = (0,1,0) direction and used imaging grid 3.

For transects, file names are of the format [map type]_[final basin diameter]_[mu]_[thresh]_[number of reversal histories].txt where mu and thresh correspondto those values for the Poisson distribution used to generate the random reversal history (see manuscript Methods). All transects assume a 200 km mapping altitude. All maps assume a magnetization intensity of 1 A/m.

For Maps, file names are of the format [map type]_[final basin diameter]_[mu]_[thresh]_[number of reversal histories]_[imaging height (km)].txt 


Possible imaging grids are defined as:
ig3=np.meshgrid(np.arange(-40e3,41e3,10e3),np.arange(-1250e3,1251e3,20e3))
ig5=np.meshgrid(np.arange(-750e3,751e3,25e3),np.arange(-750e3,751e3,25e3))
ig6=np.meshgrid(np.arange(-500e3,501e3,15e3),np.arange(-500e3,501e3,15e3))
ig9=np.meshgrid(np.arange(-1250e3,1251e3,50e3),np.arange(-1250e3,1251e3,50e3))
